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Abstract 

Staggered chiral perturbation theory (SXPT) takes into account the "fourth-root trick" for re- 
ducing unwanted (taste) degrees of freedom with staggered quarks by multiplying the contribution 
of each sea quark loop by a factor of 1/4. In the special case of four staggered fields (four flavors, 
np = 4), I show here that certain assumptions about analyticity and phase structure imply the 
validity of this procedure for representing the rooting trick in the chiral sector. I start from the 
observation that, when the four flavors are degenerate, the fourth root simply reduces np = 4 
to np = 1. One can then treat nondegenerate quark masses by expanding around the degener- 
ate limit. With additional assumptions on decoupling, the result can be extended to the more 
interesting cases of rip = 3, 2, or 1. A apparent paradox associated with the one- flavor case is 
resolved. Coupled with some expected features of unrooted staggered quarks in the continuum 
limit, in particular the restoration of taste symmetry, SXPT then implies that the fourth-root trick 
induces no problems (for example, a violation of unitarity that persists in the continuum limit) 
in the lowest energy sector of staggered lattice QCD. It also says that the theory with staggered 
valence quarks and rooted staggered sea quarks behaves like a simple, partially-quenched theory, 
not like a "mixed" theory in which sea and valence quarks have different lattice actions. In most 
cases, the assumptions made in this paper are not only sufficient but also necessary for the validity 
of SXPT, so that a variety of possible new routes for testing this validity are opened. 

PACS numbers: 12.38.Gc, 11.30.Rd, 12.39.Fe, 11.15.Ha 
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I. INTRODUCTION 

Staggered quarks |j] employ an incomplete reduction of the lattice doubling symmetry, 
and therefore have an extra degree of freedom called "taste." In four dimensions, a single 
staggered field on the lattice produces four tastes in the continuum limit. It is possible to 
interpret taste as physical flavor (u, d, s, c) by explicitly breaking the continuum taste sym- 
metry with general mass terms |2|. However, that approach leads to a variety of problems 
including complex determinants, violations of chiral symmetry even in the limit of vanishing 
light quark masses, and the necessity of fine tuning. The current standard approach — and 
the one assumed in this paper — is to introduce a separate staggered field for each physical 
flavor, and then attempt to eliminate the unwanted taste degree of freedom by taking a root 
of the staggered fermion determinant. This procedure was proposed by Marinari, Parisi and 
Rebbi j^] in a two-dimensional context; a fourth root is required in four dimensions. Such 
"rooted" staggered quarks have been used by the MILC collaboration for recent dynami- 
ca. stations Q. which give good agreement with experiment for many simple hadronic 
quantities [5|. 

There is wide-spread agreement that, whatever their practical problems in reproducing 
the desired four-flavor mass spectrum, "unrooted" staggered fermion quarks are a consistent 
way to simulate four degenerate tastes of quarks in the continuum limit. But the correctness 
of the fourth root trick to reduce four tastes to one has not been proven, and there are 
concerns expressed in the literature about its use in lattice QCD simulations 

Baa . The 

difficulties arise from the fact that taste symmetry is broken at order a 2 , where a is the 
lattice spacing. This prevents one from implementing the rooting simply by projecting the 
local four-taste staggered Dirac operator onto a local operator in a single-taste subspace. 
Without a local Dirac operator, usual physical properties of a lattice theory such as unitarity 
or universality are called into question. In the past few years, many authors have addressed 
the issue of the validity of the fourth-root procedure |£J, llO, [ill, |l2|]. While a proof is 
still lacking, the result of these investigations is to make it rather plausible that staggered 
quarks with the fourth- root trick do in fact have the correct continuum limit. At finite 
lattice spacing, however, it seems clear that the fourth root procedure introduces a variety 
of unphysical sicknesses. This follows not only from the renormalization group approach 



introduced by Shamir but also from the staggered chiral theory, as discussed below. The 
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issue is then to prove that these unphysical effects disappear or decouple in the continuum 
limit. 

Here, I start with a simpler, but related, problem: What is the chiral theory that correctly 
describes rooted staggered quarks? Lee and Sharpe found the chiral theory that corre- 
sponds to a single unrooted staggered field. In the current terminology, this is a one-flavor 
case, with four tastes. It was generalized to more than one flavor (more than one staggered 
field) and called "staggered chiral perturbation theory" (SXPT) by Aubin and Bernard jl^ . 
Certain, rather noncontroversial, assumptions go into these derivations. In particular, one 
needs to know the Symanzik theory that describes unrooted staggered quarks as one 
approaches the continuum limit. In deriving the Symanzik theory, one assumes that the 
taste, Lorentz, and translation symmetries become exact in the continuum, and that the 
lattice symmetries fit inside the continuum group in a straightforward way. In addition to 
theoretical understanding of why this should be the case H 111 Q> ll^L there is good nu- 
merical evidence for the restoration of these symmetries . To find the chiral theory for 
rooted staggered quarks, an additional assumption is needed. In Ref. it was proposed 
that one could represent the effects of the fourth root by locating the sea quark loops in 
SXPT, and then multiplying each one by a factor of 1/4. Here, I take this prescription as 
defining what I mean by SXPT for rooted staggered quarks. The question then becomes: Is 
SXPT the correct chiral theory? 

In this paper, I show that the validity of SXPT follows from certain nontrivial assumptions 
on the phase structure and mass dependence of the theory. These assumptions will be 
introduced as needed; the most important of them are also collected in the concluding 
section. While I will try to argue from simulations and experience for the plausibility of 
these assumptions, significantly more work is required to prove and/or numerically test 
them. On the other hand, in most cases it will be clear that the assumptions are not only 
sufficient for the validity of SXPT but also necessary. Tests of the assumptions therefore 
provide new means to test SXPT itself. 

Note first of all that SXPT for rooted quarks does show unphysical effects at nonzero 
lattice spacing. In the published literature, this is seen most clearly in Prelovsek's calculation 
of the flavor nonsinglet scalar correlator On the lattice, she finds intermediate-state 

contributions with mass below that of the lightest physical intermediate state (r]7i). I call this 
sickness a "unitarity violation" at finite lattice spacing, since it is due to contributions from 
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"extra" light mesons of various tastes, which only cancel completely in the continuum limit. 1 

The flavor-singlet scalar correlator provides another example of such unitarity violation. It 

las recently been worked out for the three- and one-flavor cases, both unrooted and rooted 
"1 

191 ] . Because the one-flavor case is a key test of the ideas discussed in the current paper, 
I present some relevant details in Sec. IV11 The scalar correlator at nonzero lattice spacing 
has intermediate-state contributions from light pseudo-Goldstone pions, even though a one- 
flavor theory should have only a massive pseudoscalar, the rj' '. Nevertheless, these unphysical 
states decouple from the correlator in the continuum limit. 

Thus SXPT captures some sicknesses expected of the rooted theory at nonzero lattice 
spacing. But it is not obvious that SXPT captures all such sicknesses. Perhaps there are 
other violations of unitarity, or indeed other more subtle features of the rooted theory, that 
should be present in the corresponding chiral effective theory but are missed by SXPT. I 
argue below that no such effects are missed. The starting point is a special case in which there 
is virtually no doubt about the correct chiral theory: a rooted theory with four degenerate 
quark flavors. In this case, the fourth-roots of the four determinants are identical, so their 
product just gives the determinant of a single, unrooted staggered field. (Note that the 
staggered determinant is positive, and the algorithmic treatment of the rooting trick in 
the simulations gives the positive fourth root [20].) With the noncontroversial assumptions 
mentioned above, the corresponding chiral theory is just the SXPT of Lee and Sharpe 

One can then expand around the degenerate case to treat the case of nondegenerate 
masses. For technical reasons, this requires the use of a partially quenched chiral theory with 
valence masses degenerate with those of the sea quarks. Golterman, Sharpe and Singleton 
(GSS) H Show t ha« t he phas e Structure of a ooencKeo chiral t he 0ry can be subtle, and 
analogous questions can be raised about the partially quenched theory. The use of partial 
quenching in this paper seems to be safe from any GSS subtleties. However, since the theory 
has not yet been investigated in detail using the GSS methods, I highlight a few places where 
complications could conceivably enter. Further investigation along the lines of Ref. |2]J is 
planned. 

The completion of the argument for four nondegenerate flavors requires nontrivial as- 



1 One might be tempted to describe this sickness as kind of "nonlocality" at finite lattice spacing, because 
the correlator decays at long distances at an unphysical rate. I prefer to avoid that terminology, because 
its connection with the standard issue of the locality of a Dirac operator on the lattice is indirect. 
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sumptions about the analytic structure of the mass dependence. In particular, I need to 
assume that there is no essential singularity at zero degeneracy in the difference between 
SXPT and the putative correct chiral theory. Phase transitions in the chiral theory at 
nonzero quark mass differences would also be dangerous, although the existing simulations 



4] 



can be put forward as evidence against such phase transitions, at least in the region of 
parameter space investigated to date. 

To move to the phenomenologically more interesting case of three light flavors, the mass of 
one quark can be taken large. In SXPT, it is quite clear that the heavy quark will decouple, 
leaving three-flavor SXPT. However, in the lattice QCD of rooted staggered quarks, the 
nature of the decoupling, while in my opinion plausible, requires an additional assumption. 
With this assumption, it follows that SXPT is the correct chiral description of the rooted 
three flavor theory. The process can then be repeated, leading to statements about the two- 
and the one-flavor theories. 

If SXPT is accepted as the correct chiral description, it provides strong evidence that 
rooted staggered quarks have the desired continuum limit, in other words that they are in 
the correct universality class. The point is that SXPT automatically becomes continuum 
chiral perturbation theory (XPT) in the continuum limit, modulo the usual assumptions 
on the restoration of taste symmetry in the continuum limit of unrooted staggered quarks. 
Therefore, this line of reasoning says that the low energy (pseudoscalar meson) sector of 
lattice QCD with rooted staggered quarks is, in the continuum limit, indistinguishable from 
that of ordinary QCD. This would significantly reduce the "phase space" for any possible 
sicknesses of rooted staggered quarks in the continuum limit. 

Another consequence of the arguments in this paper is more technical: If SXPT is valid, 
the lattice theory with rooted staggered sea quarks and ordinary staggered valence quarks 
(the theory in the MILC simulations ^J) behaves like a "partially quenched" theory. 2 Effec- 
tively, this means that there are symmetries that connect valence and sea quarks. As usual 
for a partially quenched theory, such symmetries may be broken in a controlled way by mass 
differences between valence and sea quarks. However, the symmetries are not broken by 
lattice corrections. The theory therefore does not behave like a "mixed" theory, in which 
valence and sea quarks have different lattice actions. In the mixed case, there are no sym- 



2 In the continuum limit, this was anticipated in Ref. [2^ . 
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metries at finite lattice spacing that connect valence and sea quarks. The chiral descriptions 
of mixed theories 23| thus have terms — vanishing in the continuum limit — that violate 
such symmetries. These terms can, for example, lead to mass splittings between mesons 
composed of two valence quarks and those composed of one valence and one sea quark. I 
show here that the chiral theory for staggered valence and rooted staggered sea quarks does 
not have such terms; corresponding valence- valence, valence-sea, and sea-sea mesons are 
degenerate. 

The remainder of this paper is organized as follows: In Sec. |Hj I discuss the replica 
trick in SXPT; this is a systematic way to find sea quark loops in the chiral theory and 
multiply each by a factor 1/4. Section HTT1 then introduces the notation needed to describe 
the various theories considered here, at both the QCD and the chiral levels, and makes some 
introductory comments about these theories. The details of my assumptions and arguments 
for SXPT in the four-flavor case are presented in Sec. IIVI while the extension to three or 
fewer flavors is treated in Sec. Section IVII shows in some detail how the one-flavor case 
works. I resolve there the apparent paradox of light pseudo-Goldstone mesons appearing the 
one-flavor chiral theory. Consequences of my arguments for the rooted theory at the QCD 
level are described in Sec. IVHI Finally, I review the assumptions and conclusions and make 
some additional remarks and speculations in Sec. IVII II 

II. REPLICA TRICK 



In SXPT for rooted staggered quarks, one needs to identify the presence of sea quark loops 
in various meson diagrams, and multiply each such loop by a factor of 1/4. The sea quark 
loops were located in Ref. [lj] by using the quark flow approach j^J. While quark flow gives 
a rather intuitive physical picture, it suffers from the disadvantage in the current case that 
it is formulated as a series of rules for tracing flavor indices, not as an algebraic statement. 
The replica trick provides an alternative approach that is systematic and algebraic. It was 
applied to partial 
for SXPT in Ref. 



y quenched theories by Damgaard and Splittorff [25] and was first used 

The replica procedure for rooted SXPT is very simple: One starts by replicating the sea- 
quark flavors, replacing each dynamical staggered field by ur identical copies, where hr is a 
(positive) integer. One then calculates straightforwardly order by order in the corresponding 
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(unrooted) SXPT, keeping the tir dependence explicit. Finally, one sets tir = 1/4. 

Note that, at any finite order in SXPT, the tir dependence is polynomial: It just comes 
from the sum over the sea quark indices in meson loops. Therefore, the process of taking 
hr — > | is straightforward and unambiguous order by order. As always in chiral perturbation 
theory, we treat the low energy constants (LECs) as free parameters for each tir. We should 
not use any relations that hold only for special values of tir — analogous to those discussed 



by Sharpe and Van de Water 27] - to reduce the number of chiral operators. If it turns 
out that we are left with some redundant operators when tir — > 4, we can always redefine 
the LECs to absorb the redundancy at the end. Within chiral perturbation theory, we do 
not worry about (nor do we have any control over) the dependence of low energy constants 
themselves on tir. Such dependence, coming from an underlying QCD-like theory, would in 
fact be nonperturbative in the strong coupling as and probably not polynomial in ur. 

At the QCD level, it is difficult to give the replica trick any meaning beyond weak- 
coupling perturbation theory, in which the ur dependence is again polynomial. Within 
weak-coupling perturbation theory, the replica trick is in fact somewhat useful, because it 
provides a convenient way of keeping track of sea-quark loops. This can aid in clarifying the 



argument in Ref. 22j of the validity of the fourth-root procedure in perturbation theory, and 
will also be helpful in Sec. IVH Bl Nonperturbatively, however, even if we were to assume 
that the Ur — > \ limit should be taken by analytic continuation, the replica trick would be 
ambiguous since there is no unique continuation from the integers. 

A related comment is that the use of the replica trick for a chiral theory is valid, a priori, 
only for order by order calculations in chiral perturbation theory. We have no guarantee of 
its correctness in general nonperturbative chiral calculations, such as the determination of 
the correct vacuum state. However, in the degenerate four-flavor theory, we know the chiral 
theory (and hence the appropriate phase) independent of the replica trick. As I move away 
from the degenerate limit, I will in any case need to assume that dependence on quark mass 
is smooth and no phase change occurs (see Sec. I1VJ) . Thus, there is no further restriction 
coming from the perturbative nature of the replica trick. 
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III. THEORIES CONSIDERED; NOTATION 



We need some notation to refer to the various versions of QCD and their corresponding 
chiral theories. Define a version of lattice QCD by (rip, np, n^pQCDi where np is the number 
of flavors (the number of staggered fields), np is the number of tastes per field, and n R is 
the number of replicas. The corresponding chiral theories are denoted by {np, rip, n R ) x . If 
riR is trivially equal to 1 (because the replica trick is not relevant), it is omitted. 

When rip, np) x or (np, np, np^pQCD are used in equations, I am referring specifically 
to the generating functionals for these theories, with sources to be discussed below. 

I focus primarily on three versions of QCD, and four versions of chiral theories: 

• (l,4) x and (1,A) L q C d'- These are the chiral and QCD theories of a single staggered 
field (one flavor) with four tastes. By (noncontroversial) assumption, the chiral theory 
(l,4) x is just the SXPT of Lee and Sharpe [l^ . No rooting is done at the QCD level, 
and no replica trick is necessary at the chiral level. 

• {rip, 4, np) x and (np,4,nR)LQCD- These are the theories for np staggered fields (np 
flavors), each replicated n R times. When % is indicated explicitly, as in this case, 
it is taken to be an integer only; no rooting is done. The chiral theories {np, 4, np)^ 
are — again by noncontroversial assumption — just those of Aubin and Bernard [l^ 
for integer (np ■ n R ) number of flavors. They are obtained from the np flavor chiral 
theories by replicating the sea-quark degrees of freedom in the chiral fields. 

• {rip, "l") x and {np, u I")lqcd- These are the chiral and QCD theories of np staggered 



fields (uf flavors) with the vDet taken at the QCD level to reduce 4 tastes to 1 for 
each flavor. Since I do not want to assume here that the rooting procedure is correct, 
I write the 1 for tastes in quotation marks. Then "l") x is by definition the chiral 
theory generated by {np, "1")lqcd- The main point of this paper is to construct 
{np, "l") x unambiguously. 

(np,4, |) x : This is the chiral theory (n F ,A,np) x , now implementing the replica trick 
by taking n R — > t , with the goal of describing rooted staggered quarks. In the liter- 



ature (e.g., Ref. [14], M, |28j, 29, 30]), it is assumed that this procedure produces the 



right chiral theory; in other words, it is assumed that (np, "l") x = (np,A, \) x . Here, 
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I define SXPT as (jipjA, 4 ) x , and then ask the question of whether SXPT is indeed 
the correct chiral theory. Note that I avoid reference to corresponding QCD theories 
"(np,4, \)lqcd" because, as discussed in Sec. HH I do not know how to give unam- 
biguous meaning beyond perturbation theory to the replica trick for those QCD-level 
theories. 

For my arguments, the chiral theories (np, 4, Wr) x are key objects. On the other hand, the 
corresponding QCD theories (ti.f,4, 71r)lqc£>, in particular (4, 4, riR)LQCD, are introduced 
for convenience, because they allow one to keep track more easily of the factors of that 
relate valence- to sea-quark matrix elements (see Sec. II V|) . These QCD-level theories can be 
eliminated at the expense of a somewhat less intuitive argument at the chiral level, related 
to quark flow. An outline of such an alternative argument is given in Sec. IIVBI it does 
however seem to require a weak additional assumption. Because the (4,4, u^lqcd theories 
are just used formally, it is probably unnecessary that the standard, broken realization of 
chiral symmetry assumed in (4, 4, Ur) x actually occurs in (4, 4, Ur)lqcd- The unpleasant 
fact that asymptotic freedom (and presumably spontaneous chiral symmetry breaking) is 
lost for tir > 1 in (4, 4, ur)lqcd seems to be irrelevant. An easy way to see this is to realize 
that the precise correspondence between (4, 4, n R ) L Q C D and (4, 4, Ur) x can be maintained 
by an artifice, 3 as follows: Note first that the order of the polynomial dependence on hr 
is bounded at a given order in chiral perturbation theory. This means there is maximum 
value of tir, n^ ax , that need be considered in order to determine the polynomial completely. 
One can then simply imagine increasing the number of colors sufficiently to ensure that 
the QCD theory has the standard, spontaneously broken, realization of chiral symmetry for 
any tir < n^ ax . Recall that the mesonic chiral theory generated by a given (4, 4, tir) lqcd 
is independent of the number of colors as long as the phase is unchanged. The numerical 
values of the LECs do depend on the number of colors, but we are uninterested in those 
values here. 

In the next section, I argue that the replica trick produces the correct chiral theory in 
the four-flavor case. In other words, I claim that 

(4,"l") x =(4,4,i) x . (1) 

3 I thank Urs Heller for this comment. 
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This should be taken as a statement about the generating functionals of the two chiral 
theories. I use "=," rather than "=," to compare two chiral theories, because what I mean 
is that they are the same functions of the LECs: True equality would only result if we 
adjusted the LECs to be the same. 

One also needs to be careful about what sources (equivalently, external fields) one is 
allowing in the Green's functions on both sides of equations such as Eq. (JIJ. For example, 
there are more sea-quark fields available in the (4, A,ur) x theory, from which (4, 4, \) x is 
obtained, than there are in the (4, "l") x theory. Unless explicitly stated otherwise, such 
generating functionals should be taken to describe partially quenched theories, with sources 
coupled to valence fields only. Ghost (bosonic) fields, degenerate with the valence fields 
and required to cancel the valence determinant, are also implicit. When I need to make 
the sources explicit, I will include any valence sources a among the arguments, for example 
(riF,nT,n R ; <t) L qcd- Identical staggered valence fields with identical valence sources are 
always assumed on both sides of equations relating generating functionals. 

IV. DETAILS OF THE ARGUMENT FOR FOUR FLAVORS 

The key ingredient is the observation that, when there are four degenerate flavors (four 
staggered fields with equal masses), the rooting procedure clearly reduces the four flavor 
theory to a one flavor theory. In other words, instead of acting on tastes and (presumably) 
reducing the four tastes per flavor to one taste per flavor, we can think of the rooting in this 
case as acting on flavor and reducing four fields to one, without affecting the tastes. Let the 
quark mass matrix be M.. The condition of degeneracy is M. = ml, where m is a number 
and I is the identity matrix in flavor space. It then follows that: 



The last equivalence in Eq. (j3J) is manifest order by order in SXPT: Since the result for any 
physical quantity is polynomial in the number of degenerate flavors, taking 4n# degenerate 
flavors and then putting = 1/4 gives the same chiral expansion as a one-flavor theory. 

One can make a stronger statement than Eq. © by adding sources and computing spe- 
cific Green's functions in the degenerate case. In order to keep the arguments simple, I 




(2) 
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generally use only taste-singlet scalar sources, which are all that are necessary to allow us 
to move beyond the degenerate mass limit. For writing explicit terms in the chiral theory, 
however, it will be convenient below to include pseudoscalar sources temporarily, since it 
is linear combinations of scalar and pseudoscalar source that transform simply under chiral 
transformations. One can also easily generalize to sources of arbitrary taste if desired. 

I start by adding introducing the scalar sources into the sea-quark sector of the QCD- 
level theory (4, "1")lqcd- Let ^i(x) be the sea quark field of flavor i at space-time point 



x. For convenience, I work in the taste representation |17| . with taste (and spin) indices on 
\1/ implicit, but there is no reason why one cannot work directly with the one-component 
staggered fields instead. The source s(x) is taken to be a Hermitian matrix in flavor space. 
The mass and source terms are then: 

m^i{x)^i(x) + ^i(x) s ij (x) Vjix) , (4, "1") case, (4) 

where sums over flavor indices i, j are implied. 

One needs to state precisely here what is meant by a rooted staggered theory with sources. 
In this paper, I always mean: (1) introduce the sources into the corresponding unrooted 
theory; (2) integrate the sea quark fields to get a determinant that is a function of the 
sources; (3) replace the determinant by its fourth root. Derivatives with respect to the 
sources, if desired, are taken only after step (3). 

Now introduce the same sources into the replica QCD theories (4, 4, tiii)lqcd, with the 
specification that a given source couples equally to all replicas. We have: 

m%(x)%(x) + %{x) s ij (x) ^j{x) , (4, 4, n R ) case. (5) 

Sums over the replica index r = 1, 2, . . . , as well as the flavor indices i and j, are implied. 

When the sources are nonzero (which includes the case of nondegenerate quark masses 
as a special case), we do not yet know that (4,4, |) x is the right chiral theory. One could 
imagine that there are extra terms in (4, "l") x that vanish in the limit s = 0. So I define 
the difference to be an unknown functional V[s]: 

(4,"l"; S ) x = (4,4,i; S ) x + \/[ S ], (6) 

where V[0] = 0. As far as we know at this point, V[s] could be quite sick. For example, 
it could generate Euclidean correlation functions with unphysical decay rates (unphysical 
intermediate states), even in the continuum limit. 

11 



There are further restrictions on V[s] coming from the fact that the two chiral theories 
must be equivalent when there is exact flavor symmetry. We must have V[s] —0 whenever 
s(x) is proportional to the identity in flavor space or can be brought there by an SU(4)l x 
SU(4)r chiral flavor rotation. Therefore it takes some care even to write down a possible 
term in V[s}. 

I temporarily add a Hermitian pseudoscalar source p(x) to the theories. For example, 
corresponding to Eq. (jlj) is 

Mi(x) KP ij (x) #i(aO , (4, "1") case. (7) 

The spurion combinations h = ml + s + ip and h* = ml + s — ip transform simply under 
chiral rotations L G SU(4)l and R G SU(4:)r: 

h^LhR^, rf^RhtlS. (8) 

If 

h(x) = c{x)U , h)(x) = c*{x)U^ , (9) 

where U G SU(4) is a constant matrix and c(x) is a c-number function, then h(x) and h'(x) 
can be made everywhere proportional to the identity by the chiral rotation R = U, L = I 
and there is exact flavor symmetry, unbroken by masses or sources. 

We can now look for possible terms in V, at first expressed as functionals of h and . 
An example that satisfies the above requirements is 

Vi= f d 4 xd 4 y (n^W) (Tr [h(x) h\x) h(y) h\y)] - ^Tr [h{x) h\x)] Tr [h(y) h\y)] ) 

J \ / x,y 

(10) 

where Tr is a flavor trace, and 1/M a distance scale that might not go to zero in the 
continuum. For example, one could have M = kAQCD, where k is some constant. In the 
worst case, M might not even correspond to the mass of any physical particle in QCD. 

Removing the pseudoscalar source p(x) and keeping only the lowest nonvanishing term 
in s, one gets the following example of a possible contribution to V[s]: 

V 1 = Am 2 J d 4 x d 4 y (^j^J [s{x)s(y)] - |Tr \s{x)\ Tr \s{y)\ ) (11) 

The goal is of course to prove that V[s] actually vanishes. 
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A. Expansion around the Degenerate Theory 



If we take derivatives of the generating functional with respect to s and evaluate them 
at s = 0, we will have Green's functions for degenerate quark masses. At the level of the 
chiral theories, I claim that Eqs. (J2J) and (jSJ) (modulo some technical assumptions) actually 
imply the stronger statement: 











FT - . . . , (4, "1"; s) x = [T A ±A,\\ s) x (12) 



s=0 



ds in ^ n (x n ) 

for any given combination of derivatives with respect to s. 

A difficulty in proving Eq. ()12|) is that, as soon as the sources are taken to be nonzero 
in order to compute the derivatives, we no longer know that (4, "l") x and (4,4, \) x are the 
same. Further, I must avoid the use of (4, 4, jjiQCDy which is not well defined. Finally, 
I cannot use (1,4:)lqcd easily as an intermediate step, because sea quark sources with 
nontrivial flavor (s^) cannot be inserted into a one-flavor theory. 

The need for nonzero sea-quark sources in Eq. (|12J1 can be circumvented by using valence 
sectors, in other words, by considering the partially quenched version of Eqs. (j2j) and Q. 
I thus introduce into all theories of interest an arbitrary number, ny, of staggered valence 
fields q a , where a = 1, 2, . . . ny is the valence flavor index. These have degenerate mass fh 
and are coupled to valence sources a a/ 3, giving mass and source terms as follows: 

fnq a (x)q a (x) + q a (x) a af5 (x) q p {x) , (13) 

with sums over a and j3 implied. The valence-quark source a alS is exactly analogous to the 
sea-quark source they only differ in the type of quarks to which they couple. 

I also introduce ny corresponding ghost (bosonic) quarks, again with degenerate mass 
fh. These ghosts do not couple to the a alS source, so that derivatives with respect to <r af3 
produce Green's functions made purely of (fermionic) valence quarks. When a a/3 = 0, the 
valence and ghost determinants cancel. 

The partially quenched version of Eq. (|2j) remains valid, since the valence/ghost sectors 
are identical on both sides, and the sea-quark determinants are equal as long as the sea-quark 
source s vanishes (giving degenerate masses): 

(4,"l"; S = 0, ( x) LQCD = (1,4; s = 0, a) LQCD , (14) 

where sea and valence sources are indicated explicitly. 
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The equality of generating functional must also be true for the corresponding chiral 
theories: 

(4,"l"; S = 0,a) x = (1,4; s = 0,a) x , (15) 

This follows by definition of what it means to be the corresponding chiral theory. I am 
assuming that such partially quenched chiral theories exist. But note that the starting LQCD 
theories both have local actions, so this appears to be a rather safe assumption. I am not 
claiming, however, that I know explicitly how to calculate ghost or valence Green's functions 
in either of these chiral theories. My expectation is that the "naive" meson Feynman rules, 
which follow from the methods of Ref. (2^, are probably correct. However, to prove that 



would require an analysis along the lines of Ref. |2l| to determine the proper saddle point for 
the mesons constructed from valence or ghost quarks, around which the chiral perturbation 
theory can be developed. Such an analysis is in progress. 

In discussing Eq. (J3J), I claimed that the equivalence of the (l,4) x and (4,4, \) x theories 
is "manifest" order by order in SXPT. In the presence of valence/ghost fields and sources, 
the corresponding statement is almost certainly still true. Even if the saddle point for 
ghost mesons (or valence) mesons is nontrivial, it is very difficult to see how it could be 
affected, order by order, by the difference between having one sea-quark flavor or having 
degenerate sea flavors and then putting ur = 1/4. Combined with Eq. (fT5|) . this gives 

(4,"l"; S = 0,a) x = (4,4,i;s = 0,<7) x . (16) 

In the limit s = = cr, all quarks, both valence and degenerate. This means 

one can relate Green's functions constructed from sea-quark fields to those constructed from 
valence fields, or equivalently, relate derivatives with respect to s to those with respect to 
a. This is not completely straightforward, however. In the (4, "1") theory, derivatives with 
respect to s bring down factors of 1/4 from y/Det(D + fh + s) = exp \tr\n(D + fh + s). 
When more than one contraction (more than one term resulting from the derivatives) is 
possible, different contractions will be associated with different numbers of factors of 1/4. 
The power of 1/4 is just the number of quark loops implied by the contractions. 

On the other hand, with arbitrary ny, we can always adjust the flavors of the valence 
sources being differentiated so that only one contraction is possible. This means we can 
always write an arbitrary derivative of the generating functional with respect to s as a linear 
combination of derivatives with respect to a, each term being multiplied by (|) L , where L 
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is the number of valence loops in the term. The following two examples should clarify what 
I mean; take flavors i 7^ j and a 7^ (3 and do not sum over repeated indices: 



d d 
dsv{x) ds^(y) 



(4, "r;s,a = 0) LQCD 



s=0 



-- (tT[G j (x,y)G i (y,x))} 



d 



d 



4 da a P (x) daP a (y) 



—-(4, u r;s = 0,a) LQC n 



CT=0 



(17) 



d d 
ds n (x) ds ll (x) 



(4, "1"; s,a = 0) 



LQCD 



s=0 



~ (tT(Gi(x,y)Gi(y,xj)) + Qj (tr (g { (x, ar))tr (^(y, y))) 

(4, "1"; S = 0,a) L QCD 



5 



<9 



It 



4 da a P(x)daP a (y) \A J ^ Qa (^) <9^%) 



5 



5 



tr=0 

(18) 



where Gi(y,x) is the propagator of a quark of flavor i from x to y, expectation values are 
taken in the (4, "1" ) lqcd theory with A4 = ml and vanishing sources, and the traces are over 
taste and spin indices. Note that the two sides of Eq. (|T7|) or Eq. (jTSjl are just two different 
ways of expressing the expectation value of the same combination of quark propagators, so 



no subtleties of partial quenching a la Ref. 



2l| can interfere with the equality. 



With enough derivatives with respect to s, there will always be enough repeats in sea 
quark flavor indices that more than one contraction contributes. On the other hand, since 
we have an arbitrary number of valence quarks at our disposal, we can always arrange the 
valence flavors in the derivatives with respect to a so that only one contraction occurs. 

In the (4, 4, Ur)lqcd theory, equations very similar to Eqs. (fTTj) and (fTHj) hold, with the 
simple replacement \ — ► ur. The factors of ur are produced by the sum over replicas for 
each quark loop. 

For an arbitrary k th derivative of (4, u 1")lqcd or (4, 4, n R ) L Q C D with respect to s, we 
therefore can write: 

A d 



n=l v ' 

L c k 



C 



E0n 



n=l 



(4, "l"; S = 0,a) iQCD 



(7=0 



(19) 



d 



n=l 
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en lc ii 



d 



(4,4,n fl ; s = 0,ct) 



LQCD 



(7=0 



(20) 



C n=l 

where C labels a particular contraction with Lq valence quark loops, and the valence flavor 
indices a% and (3% are adjusted to make only that contraction possible. The key point in 
Eqs. (II 9J) and (f2*U|) is that the same arrangements of valence flavor indices and powers Lq 
work in both cases. 

We now pass to the chiral theory in both cases, giving: 

k 



n 



(4, "1"; 8,a = 0\ 



■ i \ L c k 

e(j) n 



n=l 



da 01 " 13 " (x r 



(4, "l"; S = 0,a) 



(7=0 



(21) 



d 



n=l 



(4,4,n fl ; s,a = 0). 
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(4,4,n fl ; s = 0,(t), 



CT=0 



(22) 



At any finite order in chiral perturbation theory, both sides of Eq. (J22)) are polynomial in 
ur. Therefore the limit — > j is well defined: 

A: 



n 



n=l 



<9 



Qglnjn 



(4,4,i; s,<7 = 0), 



Etiin 



d 



(4,4i; s = 0,«t) 



cr=0 



(23) 



The right-hand sides of Eqs. (|2~Tj) and (J23|l are now equal by Eq. (fTtjj) . On the left-hand 
sides, the valence and ghost contributions cancel completely since a = 0, so we may eliminate 
those fields. This proves Eq. (JT2j). 



B. Assumptions and Results in the Four Flavor Theory 



Equation (|12j) . together with the definition of V[s], Eq. (jUJ), imply that V[s] and all of its 
derivatives vanish at s = 0: 

k 







n y (As'"-'» (,/•„ i 

n=l x v ' 







(24) 



s=0 
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Thus terms like V\ in Eq. (fTTj) are ruled out. Indeed, if V[s] is assumed to be an analytic 
function, 4 with any number of isolated singularities, it follows that V[s] =0 everywhere. 
In other words, Eq. (|TJ is true under this assumption. Normally one expects that when a 
function is expanded in a Taylor series around some point, the expansion will have a finite 
radius of convergence, given by the location of the closest singularity. But here, every term 
in the expansion is zero, so we can continue past any purported isolated singularity, and 
thereby show that the singularity is actually absent. 

Note that (4, 4, \) x , as a limit of the replica theories when ur — > \, is only defined order 
by order in chiral perturbation theory. By definition, therefore, the vacuum state of (4, 4, \) x 
has the standard broken realization of chiral symmetry that appears in (4, 4, n^) x . We know 
this is the correct nonperturbative vacuum in the degenerate limit, because there one can 
use the chiral theory (1, 4) x , for which no replica trick is needed. Now, if V[s] really vanishes 
everywhere, then (4,4, \) x is the correct chiral theory even for nondegenerate masses, and 
the vacuum must therefore remain the standard one. Thus the assumption of analyticity 
includes the assumption that there is no phase change in (4, "l") x as a function of s. 

Of course, the assumption of analyticity of V[s] is a nontrivial one. It could go wrong in 
two ways. First of all, there may be a connected "line" of singularities, an actual "domain 
boundary" that prevents one from extending V[s] = arbitrarily far from s = 0. Of course, 
XPT or SXPT must eventually break down for large enough quark masses, so it is meaningless 
to imagine extending Eq. (P) to mass differences that put one or more masses outside the 
range of chiral perturbation theory. But here I am talking about possible singularities that 
would prevent extending V[s] = over the whole range where SXPT applies. If such a 
boundary occurred, it would probably imply a phase change: that the true ground state for 
(4, "l") x changes discontinuously from the ground state assumed by (4,4, \) x . Although I 
cannot rule out this possibility from first principles, it seems rather unlikely that a phase 
change would produce small enough discrepancies to have escaped detection in the MILC 
simulations and their comparison with SXPT predictions j^J. But the effects of a phase 
change that occurred outside the (rather wide) range of masses or lattice spacings studied 
by MILC would probably not have been noticed. In addition, since the MILC simulations 

4 At this point it is sufficient for my purposes to restrict s to a constant matrix, just giving the mass 
differences. Therefore V can be thought of as a function, not a functional, and there is no subtlety with 
concepts such as analyticity. 
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involve three flavors, the logical possibility exists that a phase change occurs with four flavors 
but disappears when the fourth quark is decoupled. On the positive side, note that (4, 4, \) x 
automatically becomes standard continuum XPT in the continuum limit (see Sec. IV11 AJ) . 
Therefore, if V[s] 7^ outside some mass region, we must at least have V[s] —>■ in the 
continuum limit to avoid the bizarre scenario in which (4, "1")lqcd is a valid four-flavor 
QCD theory in some range of quark mass differences but not outside this range. 

A second way that the analyticity assumption could go wrong would be the presence of 
essential singularities in V[s] for all values of s such that the flavor symmetry is exact. For 
example, one could imagine that V[s] oc exp(— l/Vf). Although I cannot rule them out 
at this point, such singularities seem implausible to me, since we are expanding around a 
massive theory in Euclidean space and there are thus no obvious infrared problems. Note 
that I am not assuming that (4, "l") x and (4, 4, \) x separately are analytic in s around s = 
(or any other degenerate point), only that their difference is. In Sec. IVHH I speculate on a 
possible proof of the absence of an essential singularity in V[s] at s — 0. 

The assumption that V[s] is analytic is equivalent to the assumption that the expansion 
of V[s] around s = is convergent. The reader may therefore object that this assumption 
is too strong, since we do not expect convergent weak coupling expansions in quantum field 
theories. It is therefore useful to review why we believe that usual weak coupling expansions 
are at best asymptotic. The main reason comes from the factorial growth of large orders in 
perturbation theory (jy]. In the current case, however, the large orders terms in perturbation 
theory of V[s] in s are not growing factorially — in fact they are all zero! An alternative line 
of reasoning for QED is due to Dyson [32]. He argued that the expansion in a around a = 
must be asymptotic because a < leads to an unstable vacuum and therefore cannot be 
smoothly connected to the a>0 region. In fact, this argument has been shown to be flawed 
, since it is possible to define the theory consistently for a < and to obtain it by analytic 
continuation from a > 0. In any case, however, we have no similar reason to suspect that 
the difference of the chiral theories (or either of the chiral theories itself) becomes unstable 
as soon as non-zero mass differences are turned on. 

Of course, arguing that we have no reason to expect nonanalyticity in V[s] is far from 
proving that V[s] is analytic. This remains an assumption. Note that it can be turned 
around: if V[s] is not analytic then, from Eq. (|2*lj) . V[s] 7^0, so SXPT for four flavors must 
be incorrect. 
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As mentioned in Sec. IIII1 the QCD-level theories (4, 4, urIlqcd are use d in Sec. IIV Al for 
convenience; if desired, their use can be eliminated at the expense of an additional weak 
assumption about the partially quenched chiral theory. I now sketch that argument; the 
example presented in Sec. IV II can be used as an illustration of this kind of analysis. One 
needs to derive Eq. (|22|) directly in the chiral theories. It is not hard to see how to prove 
this at the chiral level, using a technique that is basically quark-flow analysis. Since the 
(vector) flavor and replica symmetries are exact in (4, 4,nn) x , one can always follow the 
replica indices though the SXPT diagrams, starting on one source index and continuing 
until one reaches another source index (on the same or a different source). Each such loop 
corresponds exactly to a quark loop at the QCD level and produces one factor of Ur. The 
same analysis then needs to be repeated for diagrams with valence quark indices. Note 
that this argument assumes that, at the chiral level, mesons made from (fermionic) valence 
or sea quarks have identical Feynman rules, except for the counting factors coming from 
replication. The ordinary, bosonic, symmetries relating fermionic valence and sea quarks 
should guarantee this, as long as such symmetries are not spontaneously broken in the chiral 
theories. Since a rigorous analysis of the partially quenched chiral theory along the lines of 
Ref. [2j| is still lacking, this absence of symmetry breaking must be taken as an assumption 
at this point if one wants to do without the use of the QCD-level theories (4, 4, tir) lqcd- 
However, it is difficult to see how it could go wrong. 



V. EXTENSION TO FEWER THAN FOUR FLAVORS 

The most interesting cases phenomenologically are three light flavors (u, d, s), or, at 
extremely low energies, two light flavors (u, d). To extend the above argument to np < 4, 
we can start by taking one the mass of one of the four quarks large and using decoupling 
ideas Q|. Call this quark the charm quark, with mass m c . The difficult point here is that 
the relation Eq. (Q) can only be used where chiral perturbation theory is applicable, so we 
cannot just take m c — > oo on both sides of Eq. (0) and then appeal to decoupling. 

In the real world, we know that the effective coupling of XPT for the strange quark is 
roughly M 2 K /(M 2 fl) ~ 0.2 0, with f n = 131 MeV. So it is likely that XPT breaks down 
completely for quark masses that are not very much larger than the physical strange quark 
mass, mP hys . For concreteness, imagine the breakdown occurs at ~2mP hys , in other words 
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for meson masses greater than ~ 700 MeV, which is the mass of a "kaon" made with a 
strange quark of mass 2m^ hys . I want to decouple the charm quark from SXPT before this 
breakdown occurs, say at m c ~1.5mP hys . Since there is not a lot of room between this value 
of m c and mf ys , it is useful to consider first the case where m s is significantly smaller than 
mf ys . I try to argue that this rip = 3 case is correctly described by (3,4, \) x . 

With m u , rrid, and m s all small, I increase m c to m c ~ 1.5mP hys . Modulo the assumptions 
discussed in Sec. IIVB1 the relation (4, "l") x = (4,4, \) x should continue to hold for m c in 
this range. I then integrate out (decouple) the charm quark degree of freedom from the chiral 
theory (4, 4, \) x . The procedure is completely analogous to the way the strange quark is 
decoupled from the continuum SU (3) l x SU (3)# chiral theory to obtain the SU (2) L x SU (2) R 
theory j^. Since this process is perturbative, there is little doubt that what remains after 
the charm quark is decoupled will be the Nf = 3 chiral theory, (3,4, nj )-v- 5 Nevertheless, a 
check of this assumption in SXPT would be reassuring, and is planned 36]. 

Thus I expect (4, u 1")lqcd with m c ~1.5mP hys to be described at low energy by the chiral 
theory (3,4, \) x . This should remain true as m c increases further, say until m c ~ 2mP hys , 
which is nominally the largest mass for which Eq. (JIJ) applies. 

Consider what happens to (4, "1")lqcd as rn c continues to increase beyond the applica- 
bility of Eq. (JTJ). When m c gets to be of order of the cutoff, m c ~ 1/a, one expects that it 
will decouple in the usual way from the QCD-level theory, leaving (3, u 1")lqcd- The only 
effect of the charm quark should be renormalizations of the (3, "1")lqcd couplings. The 
decoupling would be virtually certain if (4, u 1")lqcd were a normal theory described by a 
local lattice action. Because of the rooting procedure, though, there may be some doubt 
as to whether decoupling actually occurs. We can avoid this concern by increasing m c still 
further, until m c ^ 1/a. At that point, m c is much larger than all eigenvalues of the Dirac 
operator D, and {\/Det(D + m c ) becomes independent of the gauge field. Therefore the 
charm quark certainly decouples from (4, "1")lqcd, leaving (3, "1")lqcd- 

I am now ready to state the main assumption of this section: As m c is increased from 
~2mP hys to m c ^> 1/a, the low energy physics of (4, u 1")lqcd is unaffected, except perhaps 



5 One should also be close enough to the continuum that the taste-splittings are relatively small, so that a 
non-Goldstone meson made from light quarks is significantly lighter than any meson with a charm quark. 
This makes the MILC "coarse" lattice, with splittings as large as ~ 450 MeV in the chiral limit, rather 
problematic; while the "fine" lattices (largest splittings ~ 250 MeV), should be acceptable. 
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by renormalizations of the LECs. Here "low energy physics" means the physics of particles 
with masses and energies <C 700 MeV. An alternative way of stating the assumption is to say 
that Eq. JTJ continues to be meaningful as m c is increased from ~2mP hys to m c ^> 1/a, as 
long as (4, 4, |) x is interpreted to mean the chiral theory with the charm quark decoupled, 
and these theories are only used at low energy. 

I believe the assumption is plausible because the chiral theory shows that m c is already 
decoupled from the low energy physics by m c ~ 1.5m^ hys . I am simply assuming that it stays 
decoupled as its mass is increased further. 

The conclusion then follows immediately: (3,4, \) x is the correct chiral theory for 
(4, u V')lqcd a ^ m c ~ 2mP hys . By assumption, it remains the correct theory as m c is in- 
creased to ^> 1/a, at which point (4, u 1")lqcd becomes (3, u 1")lqcd- Thus 

(3,"l") x = (3,4,|) x . (25) 

Note that my decoupling assumption is not only sufficient for Eq. (J25|) . but also neces- 
sary. Any new physical effects entering in the region 2m phys <m c <l/a are automatically 
violations of the chiral theory (3,4, \) x . 

For the moment, Eq. (|23j) is only true for the three masses m u , rrid, m s <C mP hys , because 
these masses needed to be kept small in order to provide a clean decoupling when m c ~ 
1.5mP hys . A line of reasoning parallel to that in Sec lIVBl can now be applied: Once Eq. ()25|) 
is known to be valid for some range of quark masses, then the difference between the two 
theories must vanish everywhere if it is analytic. The analyticity could be violated by a 
phase boundary at some values of the quark mass differences. However, I can again point 
to the MILC simulations jj| as evidence against a phase boundary within the region of 
parameter space that has been studied. 

The arguments (and assumptions) of this section may now be repeated to show (2, "l") x = 
(2,4,|) x and (1, "l") x = (1, 4, ±) x . 

VI. RESOLUTION OF A PARADOX IN THE ONE-FLAVOR THEORY 

An interesting paradox arises from the final result of the previous section for rip — 1. 
Because of the anomaly, a theory with a single quark flavor should have no light (pseudo) 
Goldstone bosons, but only a heavy pseudoscalar, the i]'. On the other hand, the SXPT 
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for a single rooted staggered flavor contains 16 pseudoscalars ("pions"), of which only one, 
the taste-singlet, is heavy. The weightings of the contributions of the different pions in this 
rooted theory have factors of 1/4 compared to those in the unrooted, four-taste theory, but 
all the pions certainly contribute in both rooted and unrooted cases at finite lattice spacing. 
Even in pure-glue correlation functions, the light pions will appear as intermediate states. 6 
If my previous arguments are correct, then we know the chiral theory for a single flavor 
of rooted staggered quarks, and it will produce the correct chiral theory in the continuum 
limit for QCD with a single flavor. The only way out of the paradox is therefore that the 
light pions decouple from physical correlation functions in the continuum limit. In this 
section, I present a particular example to show in detail how the decoupling takes place at 
leading order in the chiral theory. This is a special case of the calculations of scalar, isoscalar 
correlators for various numbers of flavors worked out in Ref. [3], and will be discussed in 
more detail there. 

Gluon or glueball interpolating fields at physical momenta (<C 1/a) couple only to taste- 
invariant combinations of the quark fields. To mock up a pure-glue correlation function, we 
add a taste-singlet scalar source to the rooted one-flavor theory: 

Source = S(*)*(Z)*(Z) . (26) 

Here s(z) is basically the same as the sources considered previously in Eq. (J3J), except that 
there are no flavor indices since np = 1. 

The generating functional of this theory, (1, "1")lqcd, is obtained by computing the 
fermion determinant in the presence of the source s, taking its fourth root, and then inte- 
grating over gauge fields. In order to show explicitly the factors resulting from the rooting, I 
will take the R th power of the determinant, and only set R = 1/4 at the end. The generating 
functional is thus given by: 

n uv ,, I VA ex V {-S G (A) + gtr (In (D + m + s))} 

1 ' 1 ,LQCD JVA exp{-S G (A) + Rtr (In (D + m))} ' 1 ] 

where D is the Dirac operator for the staggered field, m is its mass, A represents the gauge 

fields, with action Sq(A), and DA is the gauge measure. As usual, one should imagine 

that additional valence quark fields (and the corresponding commuting ghost quark fields 

to cancel the valence determinant |37[) are included as needed. 



6 I thank Andreas Kronfeld for emphasizing to me the importance of addressing this paradox. 
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Note that R in Eq. (|27|) is a parameter appearing in the QCD-level generating function. It 
is logically independent from n R , which is the number of sea-quark replicas and is introduced 
(later) as a way of representing the rooting trick at the chiral level. Of course, in the end 
we want to set both R and n R to 1/4. 

I want to calculate 

G(x-y) = I (1, "l'jiQCD — I a / \ (1, "l'jiQCZ) I 7r~7 r (1, "1")lQCD 

V9s(x) 9s(y) v y s=0 V^s(x) v J s= o\ds(y) J s =o 

(28) 

The second term subtracts off the limit at infinite separation, proportional to (ty^) 2 . What 
remains is just the connected part of the correlation function at the QCD level. We are 
interested in the lightest particles that appear as intermediate states in the decay of G(x—y) 
at large \x — y\. 

There are two possible contractions contributing to the first term in G(x—y), as in Eq. (|18j): 
while there is only one contraction in each of the factors in the second term. Introducing 
valence quarks q a , degenerate with the sea quarks, and corresponding valence sources a a/3 , 
I rewrite the contributions in terms of valence quark contractions. With a ^ (3, one has 

" R2 " 1 " )l0CD L (s^te)' 1 ' " r)LQCD l, (29) 

Here and below the sea quark source s has been set to zero. 

The contractions in Eq. (|29|) are shown in Fig. ^ The first term (multiplied by R) is 
pictured in Fig. da); the second term (multiplied by R 2 ), in Fig. dh). Arbitrary numbers 
of gluon propagators and sea quark loops are implied, except that the third term in Eq. (J2HJ) 
is taken into account by omitting disconnected contributions to Fig. db). 

By the arguments of this paper, we should be able to calculate the low-mass contributions 
to G(x—y) using the appropriate SXPT. That theory is (1,4, tir) x , with n R set to 1/4 after 
the calculation to implement the replica trick. I append to (1, 4, ur) x the valence degrees of 
freedom associated with the two flavors a and (3 in Eq. (|29|). as well as the corresponding 
two commuting ghost quarks. Including taste degrees of freedom, the symmetry group of 



(l,4,n R ) x is then S'f/(4n i? + 8|8) L x SU(4n R + 8\8) R . Following the notation of [14|, I define 
a meson field $, which is a (An R + 16) x (An R + 16) Hermitian matrix, and the unitary 
matrix S = exp(«$//), where / is the LO pion decay constant. With a and b flavor indices, 
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FIG. 1: Valence quark contractions in the scalar propagator G(x — y), corresponding to Eq. (|29[) . 
The solid dots represent the source a. Only the valence quark lines are shown; completely discon- 
nected contributions to (b) should be omitted. 

running over both valence and sea flavor, we may write 

16 

® ab = ^ ts , (30) 

where the <3>^ correspond to mesons of specific taste and flavor, and t s are the 16 taste 
generators 

{ts} = {I, UG" > ")> Us, 6} ■ (31) 

with £ M the 4x4 taste matrices that correspond to the Dirac gamma matrices. All quarks 
(sea and valence) are degenerate, with mass m. 

G(x — y) will be calculated at leading order (LO) in SXPT. The valence source a cou- 
ples exactly like the valence mass term, giving a contribution to the LO Euclidean chiral 
Lagrangian: 

■^source = ^rp tr(Sp T + Sp T ) , (32) 

where \i is the chiral condensate, p and r are valence flavor indices (summed over valence- 
quark, but not ghost-quark, flavors), and tr indicates a trace over taste indices only. There 
are also terms quadratic in o appearing in the next-to-leading order Lagrangian; they may 
be ignored because they contribute only to contact terms in G{x — y) to the order we are 
working. 

To convert Eq. to the chiral level, we just replace (1, u 1")lqcd with (1,4, tir) x . 
Then, using Eq. and expanding E and to second order in $, we have 

G(x-y) = Rn 2 ($~ a (x) $%{x) ®%(y) + R 2 p 2 (^ a (x) $f a (x) <$>%{y) $g(y)\ 

\ / \ / conn 

(33) 
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(a) (b) (c) 

FIG. 2: Lowest order SXPT meson diagrams coming from Eq. and corresponding to Fig. ^ 

As in Fig. ^ a solid dot is a source, a. The cross represents one or more insertions of a "hairpin" 
vertex, and hence indicates a meson propagator that is disconnected as a quark-flow diagram. 

where there are implicit sums over the taste indices S and S', as well as over the (generic) 
flavor indices a and b, but not over the valence flavor indices a ^ (3. The subscript "conn" on 
the second term means that only those meson contractions that connect the source points x 
and y should be included. This restriction arises from the cancellations due to the last term 
in Eq. (|29|). (The first term Eq. ()33|) does not need a "conn" subscript because the valence 
indices require that all contractions connect x to y.) Cancellations of the disconnected pieces 
are also responsible for the absence in Eq. (fHHj) of contributions from the "1" terms in the 
expansion of £ + 

Figure El shows the LO (one-loop) meson diagrams contributing to Eq. (j^H)). The crosses 
indicate a presence of one or more "hairpin" vertices, which can appear on flavor-neutral 
meson lines. In the quark-flow sense, propagators without hairpin vertices are connected; 
while those with at least one hairpin are disconnected. (See for example Fig. 1 in the first 
reference in fl^j.) Note however that even a hairpin diagram is connected in the QCD (or 
meson) sense, since gluons connect the two quark lines. 

In SXPT, hairpin vertices are of two types: The first is due to the anomaly and affects 
only taste-singlet meson propagators. In the notation of Ref. it has strength 4i71q/3 
for arbitrary numbers of flavors. The anomaly contribution to the mass-squared of the rf 
is proportional to rriQ, with the proportionality constant depending on the total number of 
flavors (more precisely in this case, on the number of replicas tir). The second kind of hairpin 
is due to the taste- violating operators that appear in SXPT. These hairpins affect only taste- 
vector and taste-axial-vector mesons at LO, and have strength proportional to a 2 . Due to 
the explicit factors of a 2 , the contributions of such taste- violating hairpins automatically 
vanish in the continuum limit. Since I am interested in the restoration of physical unitarity 
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in the continuum, I omit the taste-violating hairpins here; they of course are included in a 
complete LO calculation [l9^ . 

We now go to momentum space. The (quark-flow) connected propagator carrying Eu- 
clidean momentum p is 



where Mo is the tree-level mass of a taste- S meson 

M| = 2fim + a 2 A s , (35) 

with An the taste splitting. All quarks are degenerate so there is no need to specify the 
flavor in Eq. (|33j) . 

Keeping only the taste-singlet disconnected meson propagator, we have 

®K-P) *vM) =8ajb8 V s6 a j5 Bf jV I (p) , (36) 

/ disc 



j _ Ami 1 1 

[P} = ~~ (p 2 + Mf) W+Mfi ' ( } 

M\ = Mf + n R ^ . (38) 



where |1 



with 



Note that, by definition, Mj is the mass of any taste-singlet meson before including the effect 
of the anomaly hairpin. Thus all 16 of the masses listed in Eq. ([35)1 . including Mj, become 
equal in the continuum limit. The rj'j, on the other hand, is the one meson that is a flavor 
(more precisely, replica) singlet as well as a taste singlet. Its mass M n > can be found either 
by diagonalizing the complete LO mass matrix including the anomaly term, or by summing 
the geometric series of hairpin interactions. 

One could take the limit m\ — » oo in Eq. (|3*7jl to decouple the rj', since it is at least as 
heavy as particles we have integrated out of the chiral theory (e.g., the p). However, I prefer 
to leave m finite so we may see explicitly how the rf remains after all the light mesons 
cancel in the continuum limit. 

I now consider the meson contractions that contribute to Eq. ()33j1 . Although it is not 
necessary to use a quark-flow picture here, since adjustment for the rooting is automatically 
taken into account by setting ur = 1/4, quark flow gives a nice physical picture. In Fig. 01 
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(d) (e) 

FIG. 3: Quark flow diagrams corresponding to the SXPT contributions of Fig. [21 Not shown are 
two additional diagrams that are very similar to (b) and (c) but have the roles of valence quarks 
a and (3 interchanged. Diagrams (a) and (d) have no hairpin vertices and correspond to Fig. [2^a); 
diagrams (b) and (c) have one hairpin vertex and correspond to Fig. E{b); while diagram (e), with 
two hairpin vertices, corresponds to Fig. |^c). In meson lines with hairpin vertices, a summation 
of sea-quark loop insertions is implied. 

I therefore show the quark-flow diagrams that correspond to various meson contractions. 
Figure Eta), (b), and (c) come from the R term in Eq. Note that, like Fig. [Tfa) 

from which they arise, Fig. Efa), (b), and (c) have a single valence-quark loop connecting 
the sources (shown by solid dots). Similarly, Fig. Ofd) and (e) come from the R 2 term in 
Eq. (jB3J), and, like Fig. IHb), have two separate valence-quark loops. 

When a and b in Eq. (jHSJ) quark flavors i and j, 7 connected meson propagators 

7 I use Latin indices from the middle of the alphabet (i, j, . . .) for sea quark flavors (replicas, here), 
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are only possible for the term proportional to R in Eq. ()33|). and require i = j. (In the R 2 
term the flavors do not match up.) This generates Fig. OK a), which is proportional to 
due to the sum over sea quark flavors j. 

When a and h are valence quark indices, connected contractions like those in Fig. 0t a ) 
are also possible, but there is a cancellation between valence quarks and ghost quarks, as 
follows in the quark-flow picture from the fact that Fig. 0K&) has a virtual loop. 

One additional contraction with only connected meson propagators comes from the R 2 
term in Eq. (|29|) when a = (3 and b = a. In the quark flow picture, this gives diagram 
Fig. 01(d), which is constructed entirely from valence quarks and therefore generates no 
factors of ur. 

Contractions with a single disconnected meson propagator are generated only by the R 
term in Eq. ([33)1 . It gives diagram Fig. 0(b) or the a (3 version when a = b = a or 
a = b = (3, respectively. Similarly, it gives diagram Fig. Ofc) or the a ^ f3 version when 
a = a, b = (3 or a = (3, b = a, respectively. These four terms, which correspond at the 
meson level to Fig. Efb), can be seen to have the same numerical value after adjusting the 
loop momentum assignment. 

Finally, the R 2 term in Eq. (}3*3"j) gives diagram Fig. Ofe) when a = a and b = (3. There is 
an overall symmetry factor of 2 in this case. 

We can now add the various contributions to G(q), the Fourier transform of G(x — y), 
resulting in: 

G( q ) = ^ 2 [ I (Rn R + R 2 ) V - 2 - 

W J (2tt) 4 | V R >^(p2 + M 2 )((p + qf + Mi) 

2R(4m 2 /3) f 1 1 

~ (p 2 + M 2 ) (( p + q ) 2 + M 2 ) \p 2 + M 2 , + (p + gf + M 2 , 

+ 2R 2 (4mg/3) 2 

(p 2 + M 2 ) ((p + qf + M 2 ) (p 2 + Afj) ((p + qf + M 2 ^ 

The first line in Eq. (139)) comes from Fig. 0^ a ) and (d), which give the Rur and the R 2 
terms, respectively; the second line, from Fig. Ofb), (c), and their a <-> (3 versions; the last 
line, from Fig. Ote). Note that the negative sign of the anomaly hairpin, Eq. ([37)1 . makes 

Greek indices (a, (3, . . . ) for fermionic valence quark flavors, and Latin indices from the beginning of the 
alphabet (a, b, ... ) for generic valence or sea flavors. 
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the second line negative and leads to the possibility of cancellations among the various light 
pions. It is not an accident that the hairpin is negative: It is required in order to give a 
positive mass to the rf when the geometric series of insertions is summed. 
Using, from Eq. (|38|). 

\p 2 + Mj) - (p 2 + Mf)j = J- [((p + g) 2 + M 2 ,) - ((p + g) 2 + M 2 )" , 

(40) 

one can rewrite Eq. (|3*9*j) in a form that shows more clearly how the continuum limit works: 

G(q) = /i 2 / - — ^— < —rr-p x— 7 x- + 

J W [n 2 R (> + M 2 J ((p + qf + M*,) 

(Rn R 2 )Y^ 1 ( AR 2R2 \ 1 

| . 2^ 2i? 2 \ / 1 + 1 

+ " ^ (p 2 + Mf) ((p + g) 2 + M 2 , ) + (p2 + M 2j ((p + g) 2 + M 2) 




Setting i? = 1/4 = n^, the last line of Eq. (j4*Tj) vanishes immediately. The light pions 
then contribute only in the second line. In the continuum limit all 16 of the light masses 
Ms become degenerate, and the two terms on the second line cancel also. The remainder, 
the first line, comes from the exchange of two heavy singlet mesons (rjj), and indeed has 
the same normalization as would be found for this correlation function using a continuum 
one-flavor chiral theory. This resolves the apparent one-flavor paradox, showing that it does 
not provide a counterexample to the arguments of this paper. 



VII. CONSEQUENCES 

A. Health of the Rooted Theory 

With the usual assumption that taste symmetry is restored in the continuum limit for 
unrooted staggered quarks, (rip, 4, n R ) x becomes ordinary chiral perturbation theory for Auf- 
n R "flavors" in the continuum limit. This follows immediately from the fact that, for a given 
combination of quark flavors, all 16 taste pions become degenerate in the continuum limit 
(before the effects of the anomaly are included, which affects only the taste singlet, flavor 
singlet meson, as always). Then taking n R — > \ order by order necessarily produces standard, 
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continuum XPT for rip flavors. All existing SXPT calculations [13I. Il4 Il8l l28l l29l l3Qt] have 
this expected continuum limit. 

The statement that (np,4, \) x is the correct chiral theory for {rip, u 1")lqcd (f° r n F < 4) 
therefore has important implications for the validity of the rooting procedure itself. Since 
SXPT becomes standard XPT in the continuum limit, the low energy (light pseudoscalar 
meson) sector of n^-flavor lattice QCD with rooted staggered quarks is, in the continuum 
limit, indistinguishable in its structure from that of ordinary n^-flavor QCD. There are no 
violations of unitarity, and no introduction of unphysical nonlocal scales. 

Of course, the chiral perturbation theory arguments presented in this paper do not address 
possible sickness due to rooting that would appear in sectors of the theory not described 
by XPT. Nevertheless, the extension of my arguments to at least some sectors other than 
that of the light pseudoscalar mesons seems possible. In particular, heavy-light physics can 
be described by the addition of a valence heavy quark with a nonstaggered action to the 
existing SXPT framework [28]. As such, the arguments in Sec. IIVI should also apply in 
the heavy-light case with rip = 4 sea quarks, implying that it too is free from unphysical 
effects in the continuum limit. Further, I see no obvious problems with an extension to 
rip < 4, since the heavy-quark can be treated by heavy quark effective theory at both the 
QCD and the chiral level, and thus does not introduce a new scale that would interfere with 
decoupling. The case of baryons, described by staggered heavy-baryon chiral perturbation 
theory |3flj . also seems straightforward for rip = 4. However, the artifice of increasing the 
number of colors at the QCD level is not applicable in this case, because it changes the 
nature of the baryons. Therefore, any counting arguments like those in Sec. IIV Al would 
need to be performed at the chiral level only. In addition, it is not obvious that one can use 
decoupling to analyze the rip < 4 cases, since we would now have the baryon mass scale at 
the QCD level between 700 MeV and 1/a. 

One caveat should be added to the discussion of this section: Since the chiral expansion 
expresses physical quantities in terms of unknown LECs, the statement that SXPT is valid 
does not by itself imply that the LECs generated by the rooted staggered theory take on their 
correct (real QCD) values in the continuum limit. On the other hand, in the four-flavor case 
we do know that the LECs are correct in the degenerate case. This follows from universality 
since the degenerate action is local. But the LECs are by definition mass independent, so if 
SXPT is indeed the right chiral theory for four nondegenerate flavors, the LECs are per force 
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also correct. With fewer than four flavors, though, my assumptions on decoupling do not 
appear to be strong enough to continue to guarantee correct LECs. For that one would need 
to show universality at the lattice QCD level (see for example Ref. or to argue from the 
agreement of rooted staggered simulations with experiment j^|. Of course, numerical checks 
against experiment are not proofs, and they run the risk, at least in principle, of confounding 
small violations of universality due to rooted staggered quarks with small violations of the 
Standard Model. Such checks will become more convincing when one can see agreement 
between at least two different lattice fermion approaches. 

B. A "Mixed" Theory? 

n 

In current dynamical staggered simulations |4| , the fourth-root trick is applied to the sea 
quarks, while the valence quarks are described by ordinary staggered fields. In this section, 
I call this situation a "rooted-staggered theory" for simplicity. Because valence and sea 
quarks are treated differently, it has been suggested Q that rooted-staggered theories fall 
into the class called "mixed," where the valence and sea quarks have fundamentally different 
lattice actions. In mixed theories the mass renormalizations of sea and valence quarks are 
different, meaning in particular that there is no simple way to ensure that sea and valence 
quarks have the same physical mass. Further, the continuum symmetries that would rotate 
valence and sea quark into each other are violated by discretization effects. This implies, 
for example, that even if the quark masses are adjusted to make the mass of a meson with 
two valence quarks equal to the mass of a meson with two sea quarks, the mass of a meson 
with one valence and one sea quark will be different by O(a) or 0(a 2 ) terms. Such terms 
show up in new operators in the XPT for the mixed theory 23]. 

I claim, however, that the rooted-staggered case is not a mixed case, but in fact resembles 
much more closely a partially quenched theory, where the symmetries between valence and 
sea quarks are violated only by explicit differences in quark masses. 

First of all, I sketch a proof that the renormalization of sea and valence quark masses are 
the same to all orders in (weak-coupling) perturbation theory. Imagine we have determined 
the mass counterterm for a valence quark up to an including a given order in perturbation 
theory. I need to show that the same mass counterterm will work to renormalize the mass on 
a sea quark line that appears as a loop inside some other diagram. Go inside the diagram, 
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and draw a box around a self-energy insertion on a sea quark line. As remarked in Sec. |HJ 
the replica trick shows that the rooting procedure simply multiplies each sea quark loop 
by 1/4 in perturbation theory, so the self-energy insertion as well as the associated mass 
counterterms on that line are all multiplied by the same overall factor of 1/4, compared to 
the corresponding self-energy insertion and mass counterterms on a valence line. Thus the 
same counterterms work in both cases. Of course there may be additional factors of 1/4 for 
any sea quark loops that appear in sub-sub-diagrams. But these will be the same for a sea 
quark line as for a valence line. 

The argument in the proceeding paragraph is based on weak-coupling perturbation theory. 
Could there be "mixed effects" that show up only nonperturbatively? I cannot answer that 
question for general nonperturbative effects, but I can answer it — modulo the assumptions 
in Sees. IIVI and |V] — for the large class of effects described by the chiral theories. The 
appropriate chiral theory is {rip, A, |) x , which is obtained order by order from {uf, 4, Wr) x . 
The latter theories have symmetries interchanging valence and sea quarks. For ny flavors 
of valence staggered quarks, the full symmetry group is in fact SU{AriRnF + 4ny|4ny)z, x 
S77(4nRnp + 4ny \Ahv)r- The taste symmetries are broken on the lattice at 0(a 2 ), but the 
"flavor subgroup" 8 U{n R nF+n v \n v )e x U{n R n F +n v \ny) r is exact up to quark mass terms. 
Extra chiral operators that would split valence-sea mesons from sea-sea or valence-valence 
mesons are forbidden by these symmetries. Since such operators are absent for all n R , they 
can have no effect when we take n R — > j. In particular, corresponding sea-sea, valence- 
valence, and valence-sea mesons are all degenerate (when the quark masses are degenerate) 
in {np., 4, 77, j r) x , and therefore in (rip, 4, \) x . Thus, at least within the context of chiral 
perturbation theory, the rooted-staggered theory behaves like a partially quenched theory, 
not like a mixed theory. 

One does have to be careful in defining the word "corresponding" in the previous para- 
graph. The valence sector of a rooted-staggered theory is "richer" than the sea sector, in that 
it includes particles in the continuum limit whose sea-sector analogues have decoupled from 
the physical theory. This is not surprising, since the purpose of the fourth root is to reduce 
four sea quark tastes to one, and there is no fourth root taken in the valence sector. A simple 

8 This flavor subgroup is described in the first paper in Ref. but is there called the "residual chiral 
group." It has been generalized here to take into account the partially quenched context. 



32 



example of this behavior can be seen from the result in Sec. IVI1 If one adds together the 
valence contractions in Eq. without the extra factor of R relating the last two terms to 
the first, then one gets a valence Green's function with no sea-quark analogue. Intermediate 
light (pseudo-Goldstone) mesons will appear as intermediate states of this Green's function 
in the continuum limit. In this sense, the rooted-staggered theory, is inherently "partially 
quenched," even in limit of equal valence and sea masses. In a normal partially quenched 
theory, one can always take more valence quarks than there are sea quarks, so one has the 
option of creating valence states that have no analogue in the sea-quark sector. The main 
difference here is that one has no choice in the matter: The physical sea-quark subspace is 
always a proper subspace of the complete theory in the continuum limit. 



VIII. CONCLUSIONS AND DISCUSSION 



Under certain assumptions that I repeat below, I have shown in this paper that staggered 
chiral perturbation theory (SXPT) correctly describes the low energy physics of four or fewer 
flavors of rooted staggered quarks. The SXPT theory (n^,4, j) x takes into account the 
fourth- root procedure by the replica trick (or equivalently, by quark- flow analysis). At finite 
lattice spacing, SXPT reproduces unphysical features of the rooting that may perhaps best 
be described as violations of unitarity, with unwanted intermediate states contributing to 
amplitudes. This is clearly seen in Ref . or in the example presented in Sec. IVII . 

Because SXPT becomes standard XPT in the continuum limit, the unitarity violations 
seen in SXPT at nonzero a must go away when a — > 0. If SXPT is indeed the correct chiral 
theory for rooted staggered quarks, then this implies that the low energy (pseudoscalar 
meson) sector of lattice QCD with rooted staggered quarks is, in the continuum limit, 
indistinguishable in its structure from that of ordinary QCD. There are no violations of 
unitarity, and no introduction of unphysical nonlocal scales. This would not, by itself, show 
that the rooting procedure is valid, because there could be problems in sectors of the theory 
not described by chiral perturbation theory. Nevertheless, it would significantly reduce the 
possible ways in which rooted staggered quarks could go wrong. 

My SXPT results also give support to the statement that the theory with staggered va- 
lence quarks and rooted staggered sea quarks is not a "mixed" theory. Like a partially 
quenched theory with the same action for the valence and sea quarks, the rooted staggered 
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theory has flavor symmetries rotating sea and valence quarks into each other. These sym- 
metries may be broken in the usual way by mass terms, but they are not broken by lattice 
corrections. 

The starting point of my argument was the observation that four flavors of degenerate 
staggered quarks simply reduce to a single flavor when the fourth root of the determinant is 
taken. To make use of this observation, I needed several assumptions, the most important 
of which are: 

1) The taste symmetry is restored in the continuum limit of the normal, unrooted stag- 
gered theory. 

2) The difference V[s] between the SXPT theory for four flavors (4, 4, \) x and the true 
chiral theory (4, "l") x is analytic in s (for space-time independent s), up to possible 
isolated singularities. 

3) As a single quark mass ("charm") is increased beyond the point at which it has de- 
coupled from the chiral theory to a scale much larger than the lattice cutoff, the low 
energy physics of (4, u 1")lqcd is unaffected, except perhaps by renormalizations of 
the LECs. 

I consider assumption 1) to be noncontroversial, and there is a lot of numerical evidence 
or it, but it has not been rigorously proven. The renormalization group methods of Shamir 



111 ] seem to me the best way to make progress on this issue. 

Assumption 2) is used in Sec. IIVI to move from degenerate to nondegenerate masses in 
the four-flavor case. The most important "obstruction" here would seem to be the possible 
existence an essential singularity in V[s] at s = 0; I speculate below on how this possibility 
might be eliminated. Note that the existence of such singularity immediately would imply 
that SXPT is incorrect. A second way the assumption could be violated would be the 
presence of a phase boundary at a finite distance from s — 0. This would imply the existence 
of a region of mass differences in which SXPT is valid, and another region of larger mass 
differences in which SXPT is invalid. Generically, I would expect an abrupt change like 
this to cause significant effects that would likely have been noticed in simulations if they 
occurred within the parameter ranges studied. Both types of potential analyticity violations 
certainly merit further investigation, however. 

34 



Assumption 3) allows me in Sec. to extend the result in the four-flavor case to the 
more interesting cases with fewer than four light flavors. It should be possible to test this 
assumption numerically by simulating a four-flavor theory in appropriate mass range and 
seeing if it is describable at low energy by the proper chiral theory with a decoupled charm 
quark, (3,4, \) x . Such tests are under consideration by the MILC Collaboration and may 
be performed in the near future. The main uncertainty is the precision at which these tests 
could be made, which would strongly influence how convincing they would be. 

To investigate a possible essential singularity, I restrict myself to diagonal sources, con- 
stant in space-time. In other words, I consider a function V of the four mass differences 
from the degenerate mass fa. To correspond with the previous notation, I write V = V(s), 
with s y = 5^€j and Sj = raj — rh. Considering s to be complex, the arguments in Sec. IIV Al 
still go through formally, although one may want to put the system in finite volume to avoid 
any dangers from j d A x with constant sources. We then have a complex function V(s), all 
of whose complex derivatives vanish at s = 0. This would forbid essential singularities in 
V(s), which do not have well-defined complex derivatives. 

What would be needed to make such an argument reasonably rigorous? On the SXPT 
side, we are defining (4, 4, \) x in (chiral) perturbation theory, so I do not expect problems at 
any finite order in adding small, complex 6j to the masses in Euclidean space. On the other 
hand, we do not know what (4, "l") x is a priori, so we would need to add ej to the masses 
in the QCD-level theory (4, u 1")lqcd- The main problem there seems to be that a complex 
s makes the determinant complex. The issue of how one chooses the phase of the fourth 
root thus becomes relevant, as it is for the case of nonzero chemical potential 38]. Unlike 
the chemical potential case, however, the imaginary part of 6j adds a constant amount to 
all eigenvalues of flavor j. Furthermore, €j may be taken very small, i.e., much less than 
both rh and Aqcd- I am hopeful therefore that any phase ambiguities can be shown to 
be manageable, but that remains to be seen. A further difficulty could come in trying to 
"match" (4, "1")lqod on t° (4,4, \) x in order to turn statements about smoothness of each 
theory separately into statements about V\s\. 

I conclude with a two additional comments: 

• Because the "wrong" mesons, which may be lighter than the physical states, contribute 
to correlation functions at nonzero lattice spacings in SXPT, the infinite-distance limit 
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of some quantities may not commute with the continuum limit. This order of limits 
issue is very similar to that concerning the chiral and continuum limits, described in 



Ref. 39]. It is not a practical problem, since of course only finite distances are relevant 
to simulations, and the extrapolation can be taken in the proper order with the aid of 
SXPT. 

There is nothing in my argument that (4,4, \) x is the correct chiral theory for four 
flavors of fourth-rooted staggered quarks that is really dependent on the fact that 
we are taking the fourth root. The same arguments would also imply that (3,4, \) x 
is the chiral theory for three flavors of staggered quarks for which the third root is 
taken! The decoupling arguments in Sec. |V| (which presumably still apply) would 
say further that (rip, A, |) x gives the chiral theory for rip < 3 flavors of third-rooted 
staggered fields. There is no contradiction here: (np, 4, g) x describes a sick theory, 
even in the continuum limit, except for the uninteresting case of 3 degenerate flavors. 
Since a staggered field always has four tastes, only a fourth (or square) root the 
root can describe an integer number of flavors (and therefore a local action) in the 
continuum limit. The np — 1 example from Sec. IVII provides a simple illustration: 
the contributions from light pions in the second and third lines of Eq. (14 lj) vanish in 
the continuum if and only if tir = R = 1/4. (I ignore the trivial case R = 0, as well 
as ur = R = — 1/4, which violates the spin-statistics theorem.) Even ur = R = 1/2 
leaves some light-pion contributions, as it should since that is really a two-flavor case. 
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